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Abstract 

The study of chaos in relativistic systems has been hampered by the observer dependence of Lyapunov 
exponents (LEs) and of conditions, such as orbit boundedness, invoked in the interpretation of LEs as 
indicators of chaos. Here we establish a general framework that overcomes both difficulties and apply the 
resulting approach to address three fundamental questions: how LEs transform under Lorentz and Rindler 
transformations and under transformations to uniformly rotating frames. The answers to the first and third 
questions show that inertial and uniformly rotating observers agree on a characterization of chaos based on 
LEs. The second question, on the other hand, is an ill-posed problem due to the event horizons inherent to 
uniformly accelerated observers. 
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The quest for an observer-independent characterization of chaos in relativistic systems yj] has 
been an intense area of research and promises to provide significant new insights into the properties 
of chaotic dynamics yj]. An important recent result [3] concerns the transformation of Lyapunov 
exponents (LEs) under spacetime diffeomorphisms. We recall that the dynamics of a bounded 
solution X(£) of a dynamical system 

is chaotic if it presents sensitive dependence on initial conditions [4]. The associated LEs Q5[] 
are given by Aj = lim sup^^ | log | 1 1, where ip^t) are solutions of the linearized equation 
4<^ = [DxF(X.(t))} ifi. Positive LEs are related to exponential divergence of initially close 
trajectories and, consequently, to chaotic dynamics. For space diffeomorphisms X = ^(Y), the 
invariance of the LEs is well established under rather general conditions (see, e.g., [f], 7]). In 
contrast, for well-behaved spacetime diffeomorphisms involving time changes of the form dr = 
A(X)dt, it has been shown [3D that the LEs transform according to 

K = Af/<A> t , (2) 

where < (A) f < oo is the time average of A along the corresponding trajectory. Therefore, 
although the values of the LEs are themselves non-invariant, their signs are preserved and assure 
an invariant criterion for chaos under spacetime transformations. This result was obtained under 
conditions for which LEs are known to be valid quantifiers of chaos, of which the most limiting 
ones are the assumptions that the system has a natural invariant probability measure and the orbits 
are bounded both before and after the transformation. 

In this Letter, we extend this result to an important class of transformations that do not preserve 
the boundedness of the orbits, and address fundamental questions on relativistic chaotic dynamics 
that require explicit in-depth investigation due to their outstanding physical properties and the 
violation of conditions invoked in the derivation of Eq. ©. The first question is how the LEs 
transform under Lorentz transformations. This question determines whether all inertial observers 
agree on a LE-based characterization of chaos. We show that the answer is affirmative despite the 
fact that the dynamics becomes unbounded with respect to at least one of the reference frames. 
We use this example to establish an extended boundedness condition for the definition of the 
LEs as indicators of chaos, which is formulated relative to the trajectories themselves rather than 
a fixed point of the phase space. The second question is how the LEs behave under Rindler 
transformations, a question equivalent to ask whether uniformly accelerated observers agree on an 



inertial characterization of chaos based on LEs. We show that this question is ill-posed because 
uniformly accelerated observers do not have access to the late-time dynamics. The latter relates to 
the fact that chaos and LEs are asymptotic concepts [8] whose definitions involve a limit t — > oo. 
We also consider transformations to uniformly rotating frames, and show that the positivity of the 
LEs remains invariant under such transformations. 

Our principal result stems from this analysis and can be stated for any system and any spacetime 
diffeomorphic transformation, as follows. For the system written in autonomous form, the LEs 
transform according to Eq. © and remain invariant indicators of chaos if, as shown below, (i) our 
extended boundedness condition is satisfied, (//) the Jacobian of the transformation is bounded, 
and (Hi) A is positive for all t and < (A) t < oo. These conditions depend not only on the 
transformation properties of the dynamical variables X and the change of reference frames but also 
on the choice of spacetime coordinates. They are automatically satisfied for global nonsingular 
transformations of bounded orbits for which inf A ±:L > whether the system is conservative, 
dissipative, mechanical, chemical, thermodynamical, electromagnetic, or fluid dynamical. These 



conditions clarify previous results [9] that seem to challenge the invariance of chaos for relativistic 
observers, and show that LEs lead to invariant conclusions about chaos. 

We first note that under a space diffeomorphism X = *(Y), system (OQ) is mapped into = 
[Dy$ (Y)] -1 F(^f(Y)), rendering the solutions of the new linearized dynamics to be related to 
those of (0Q> as cp^t) = [D Y ^(Y(t))] fl7J]. Hence, the corresponding LEs satisfy 



*-oo t \\<Pi(t)\\ t-oo t \\<Pi{t)\\ 

Suppose the solutions X(£) are limited to a compact subset of the space. Since the diffemorphism 
maps bounded solutions X(t) into bounded solutions Y(t), the matrix D Y ^ (Y (t)) is nonsingular 
and, besides, there are time-independent finite nonzero constants = sup || [Dy* (Y (t))^ 1 \ \ 
leading to 

lim - log — < Xi - Ai < lim - log L + , (4) 

t— >oo t L t— >oo t 

which imply A, = Aj [7]. This argument explores the boundedness of X(i) and Y(t) to ensure the 
existence of the constants L ± . Below we extend © and establish Eq. © for an important class of 
unbounded orbits. 

We now consider transformations of reference frame in which (OQ) describes a bounded au- 
tonomous system with respect to the initial (inertial) observers. More general transformations 
can be obtained by a composition of such transformations. We start with single-particle systems. 



While general relativity allows arbitrary spacetime coordinates, and conditions (i-iit) can be ap- 
plied to any of them, we will assume that the dynamics is described in terms of physical times 
(i.e., the time measured by observers at rest in the reference frame at the corresponding space 
coordinates). 

Lorentz transformations. We first focus on the case in which function F depends only on 
the configuration- space coordinates, such as in the evolution of a fluid element determined by a 
stream function, and consider a Lorentz boost with velocity v along the ^-direction, (ct, x, y, z) — > 
(ct 1 \x' \y' \z') = \E r " 1 (ci, x, y, z), where 

^ l (ct,x,y,z) = (-/(ct - ra/c),7(x - vt),y,z) (5) 

for 7 = 1/ a/1 — iy/c) 2 . We focus on the space spanned by the coordinates (ct, x, y, z) = (ct, x), 
where we have enlarged the configuration space in order to incorporate ct as a new coordinate. 
The extended version of (OQ) then reads 



d I 








dl I 




H 


1 F(x) / 



(6) 



where dw/dt = d(ct)/dt. 

The main advantage of this formulation is that the transformed system remains autonomous 
and the spacetime transformation can be reduced to an ordinary space diffeomorphism; it can 
be split as T o S(ct,x), where S is a transformation (w',x') = ^f^ 1 (w,x) that preserves the 
independent variable and T is a time redefinition dt! = A(w, x)dt. (Another advantage is that the 
analysis extends immediately to F with explicit time-periodic dependence.) The solutions of © 
are unbounded along the ^-direction, but this is not a problem since the nonzero LEs of system 
© are identical to those of (IB- 
There is a caveat, however: the spatial boundedness of the solutions is not preserved under 
Lorentz transformations. A trajectory confined to a bounded space-like region (sup \ \x[t) \ \ < oo) 
of the first reference frame is seen as spatially unbounded from the other inertial reference frame. 
Similar problem is observed even for Galilean transformations, but in classical dynamics one can 
adopt a reference frame where the solutions are bounded. In relativistic dynamics such a choice 
would raise questions about the invariance of the LEs, which is precisely the object of this Letter. 

To proceed we first make the crucial observation that the study of chaos can be extended to 
this class of spatially unbounded orbits, even though the same does not hold true for unbounded 
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systems in general. Indeed, sensitive dependence on initial conditions and LEs depend exclusively 
on the relative time evolution between nearby trajectories; their dependence on the reference frame 
is limited to the definition of the spacetime coordinates used to measure the distances between the 
neighboring trajectories as they evolve over identical time intervals. Therefore, chaos can be 
properly defined and LEs can be used as indicators of chaos on an unbounded trajectory y(t) 
insofar as \\y(t) — y(t)\ \ remains uniformly upper bounded for all t and all trajectories y{t) with 
initial conditions in a neighborhood ofy(0). That is, our condition is that the evolution of a small 
ball of points will remain bounded with respect to the local observers at position y(t), regardless 
of whether it remains bounded with respect to a fixed point of the reference frame. We refer to this 
as the extended boundedness condition. Note that this condition is satisfied for y(t) interpreted 
as the extended coordinates (w'(t), x'(t)) after the transformation S whenever the original system 
(OQ) is spatially bounded. 

Having shown that LEs remain valid indicators of chaos despite the spatial unboundedness of 
the transformed orbits, we now turn to the effect of the Lorentz transformations on the LEs. For 
the transformation T, from Eq. © we have 

dt' = 7 (l - ^F x (x(t)j) dt = A{x{t))dt, (7) 
where F x (x) stands for the ^-component of F(x). For \F x (x(t))\ < c, implying inf A(x(t)) > 



in the present case, we have < (A) t = lim^oo t -& = lim^oo | f* A{x{p))dp < oo ft 1 Oil . This 



allows us to factor the LEs transformed by T o S as 

Af = \l/(A) t , (8) 

where (A) t is the contribution due to T and A- corresponds to A' transformed by S. The problem is 
thus reduced to the transformation of the LEs under the spatial transformation S. The nonsingular 
nature of © assures the existence of the constants necessary to establish the bounds in © 
because, irrespective of the spatial unboundedness, the Jacobian matrix of the transformation is 
bounded. Employing the Euclidean norm to the matrix D y ^f(y) of ©, we obtain L + = L~ = 
a/(c + \v |)/(c — \v I), leading to A* = A*. In particular, all positive LEs remain positive under this 
transformation. 

Combined with Eq. ([8]), this results in A* = A'/ (A) t , which is precisely the transformation 
© previously established for the case of bounded orbits [3]. Our result does not agree with the 
result presented in |Q] for averages over local LEs but that is because that study was restricted 



to time dilatations and length contractions, which correspond to the transformation of dynamical 
variables such as volume (or the reciprocal of a density) for the time measured at a fixed point of 
the reference frame, whereas our analysis describes single-particle dynamics for the time measured 
at the position of the particle. 

If system © involves the evolution of velocities, as expected for a particle in a 3D potential, 
the Lorentz transformation © must be extended to include the transformation of u = dx/dt 
into u' = dx'/dt', which is given by u' x = r)(u x — v), u' y = r}^~ l u y , and u' z = r]^j~ l u z , where 
7] = 1/(1 — u x v/c 2 ). The resulting transformation (w, x, u) — > (w', x', u') satisfies the extended 
boundedness condition and has constants < < oo, as long as |i>| < c and < c. 

This ensures that the LEs of systems obtained by order reduction of second-order differential 
equations, which are the most common in particle dynamics, will be transformed as in © under 
Lorentz transformations. 

Rindler transformations. With respect to an inertial reference frame, an observer with con- 
stant proper acceleration a along the x direction has a hyperbolic worldline given by 



c 2 , ar , s c 2 ar 



ct(r) = — sinh — , x(r) = — cosh — , (9) 
a c a c 



where r stands for the observer's proper time. The corresponding Rindler transformation Iil2ll is 
defined by [ct, x, y, z) — > (cr(t, x),£(t, x),y, z), with 

ct(r, £) = c\ — sinh — , x(r, ^) = c\ — cosh — , (10) 
V a c V a c 

and positive £ (see Fig. [T]). The observer on hyperbole © is in the Rindler reference frame at rest 
at £ = c 2 /2a. In contrast with the Lorentz case, the Rindler transformations are non-linear in x 
and ct. 

Focusing on the space defined by the extended configuration-space coordinates, the matrix 
D y ty(y) and its inverse for the Rindler transformation of © have unit determinants but their 
largest eigenvalues diverge as c/i/2a£ cosh ar/c for £ — > 0. Therefore, one cannot identify fi- 
nite constants L ± that could be used to compare A* and A*. This behavior can be interpreted 
in terms of our extended boundedness condition, which is not satisfied in this case because 
(cr(t, x(t)),£(t, x(t)),y(t), z(t)) diverges at the light cone and is undefined beyond it. Moreover, 
from the inverse of Eqs. (fTOl) . we have 



z 2 

dr = — 

a 



r- x(t)-tF x (x(t)) 
x(t) 2 - (ct) 2 



dt = A(ct,x(t))dt, (11) 
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FIG. 1: Accessibility to the dynamics is observer dependent. The lines of fixed £ in Rindler coordinates 
correspond to hyperbolic trajectories in the coordinates (ct, x) of inertial observers. The straight dotted lines 
are the lines of constant time r. The uniformly accelerated observers are unaware of all events occurring 
in regions II and III of the original Minkowski spacetime. They only access the dynamics of a trajectory T 
during the time the trajectory crosses region 1 11311 . If the trajectory is spatially bounded with respect to the 
original observers, as assumed for system (Q]), this corresponds to an infinite time interval At but only to a 
finite time interval At. 

where A(ct, x(t)) diverges when the original solution (ct, x(t)) crosses the light cone x 2 = c 2 t 2 . 
The same holds true for the physical time dt' = ^2a^/ c 2 dr. The average (A) t is not well defined 
and, as a result, the Rindler transformed system does not have a natural probability measure against 
which the LEs could be calculated yO. Therefore, the question of how the LEs transform under 
Rindler transformations is ill-posed. 

The real origin of the problem is the horizon structure (and its counterpart structure for t — > —t) 



inherent to uniformly accelerated observers [12]. The Rindler transformation (1101) is not a global 



spacetime diffeomorphism since it maps only one quarter of the Minkowski spacetime, as shown 
in Fig.CQ Any event located above the component of the light cone corresponding to the bisectrix 
in the first and third quadrants of Minkowski spacetime will never reach the accelerated observers. 
While singularities can be an artifact of the coordinates, event horizons are an attribute of the 
reference frame. The existence of an event horizon prevents the observers from having access to 
the asymptotic dynamics of the original system. Therefore, without having access to the complete 
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dynamics, the Rindler observers cannot formulate a criterion for chaotic behavior — based on the 



observation of individual trajectories — that is valid for the original system [13Q. It is interesting 
to notice that such a problem, related to the global structure of the spacetime, manifests itself as a 
violation of our conditions for the transformation of LEs. 

If one insists on computing the LEs from a uniformly accelerated referential frame [9], one must 
note that the late-time dynamics of the extremely dilated time r — > oo does not correspond to the 
real late-time dynamics of the original system since the interval — oo < r < oo is the mapping 
of a finite time interval At. Therefore even if one could compute A[ as seen from the accelerated 
frame, this would be, in fact, a problem different from the originally proposed one. This situation 
is analogous to the limits imposed by the cosmological singularity to the determination of chaos 
in FRW cosmologies [8] and is also predicted for Rindler transformations of any other dynamical 
system and for any choice of coordinates. 

Rotating frames. The crucial role played by the event horizon in the Rindler case can be better 
appreciated if one considers a physical situation involving a non-linear transformation that does not 



introduce event horizons. This is precisely the case of uniformly rotating reference frames 

r' = r, 9' = 9 + fit, z' = z, and cdf = [g(r) + fi 2 r 2 / g(r)]dt + [VLr 2 / g{r)}d9, where g(r) = 



\/c 2 — fi 2 r 2 , fi is a constant, and t' is the physical time in the rotating frame ltl5h . This leads to 

\ c cg{r{t)) J 

where Fg(x) = d6/dt. The transformation of the LEs of © is in this case well defined since 
the extended boundedness condition is satisfied for orbits in closed sets of the physical region 
|fi|r < c for which — Vtr 2 F 9 (x) < c 2 , where both the function A(x) and the constants V s1 are 
upper and lower bounded away from zero. The latter follows from the fact that the entries of the 
Jacobian matrix D y ^{y) and its inverse for the transformation (ct, r, 9, z) — > (cf, r', 9', z') are 
all continuous for Vtr < c. A subtlety in this calculation is that in rotating frames the differential 
dt' of the physical time is not exact and cannot be integrated globally, meaning that the Jacobian 
elements involving derivatives of ct' must be determined from cdt' in the immediate neighborhood 
of a given r. The transformation t — > t' is defined locally but it can always be extended along any 
trajectory with initial condition in that neighborhood. Therefore, the LEs transform as predicted 
by (O also for the case of rotating frames. 

Generalization and discussion. Our derivation of Eq. ([8]) also demonstrates that conditions 
(i-iii) are sufficient (and usually necessary) for the validity of © in general. Indeed, while we 
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considered specific transformations and specific classes of dynamical systems in our explicit ex- 
amples, these three conditions are precisely the checkpoints we have to verify for any system 
and any transformation. The extended boundedness condition — satisfied both before and after the 
transformation in the extended space, which includes dt as an additional coordinate — guarantees 
that the system can be kept autonomous and that LEs remain valid indicators of chaos. The con- 
dition that the Jacobian is bounded — in the sense of having positive finite constants for the 
transformation in the extended space — ensures the validity of the identity A- = A*. Finally, A and 
(A) t positive and finite — again, in the extended space — guarantees that the time transformation is 
well defined and the signs of the LEs are conserved; it also guarantees that the time transformation 
is invertible, a condition we saw violated for the Rindler transformation. 

These conditions are readily applicable to any system and any change of reference frame 
and coordinates. The latter includes the choice of the time parameter or of the observers in 
the reference frame with respect to which the time is measured. In the examples above, the 
dynamical system describes the dynamics of a single particle, the dynamical variables repre- 
sents the coordinates and possibly velocities of the particle, and the time was assumed to be 
recorded locally — each time by the observer in the reference frame that is at the point where 
the particle is. However, other choices are equally valid. For a many-particle system under 
Lorentz transformation, for example, the time could be measured, e.g., with respect to the po- 
sition of one of the particles, dt' = 7 (l — j%F Xi (x(t))) dt, with respect to the center of mass, 
dt' = 7 ^1 — j% J2i Y^ln . FxX^^))^ dt, or with respect to a fixed point, dt' = jdt. Moreover, the 
dynamical system can describe physical, chemical or biological activity whose dynamical vari- 
ables do not necessarily correspond to coordinates and velocities in the physical space. In this 
general case the system can be written as J^-Xj = Fj(Xi, . . . X n ), i — 1, . . . n, and the transforma- 
tion is locally defined as (cdt, dX 1: . . . dX n ) — > (cdt', dX[, . . . dX' n ). The latter is determined by 
the change of reference frame and spacetime coordinates, (cdt, dx) — > (cdt', dx'), and depends 
on the nature of the dynamical variables, i.e., whether they transform as scalars, vectors, tensors, 
or in a different way. The choice of observers in the new reference frame is always accounted for 
through the choice of j- t x in the transformation formula dt' = (j^t'(ct, x) + V x t'(ct, x) ■ J^as) dt, 
where this term vanishes only if the time is measured (remotely) by a fixed observer. 

The results presented in this Letter address all these cases and show that, if conditions (i-iii) are 
verified, the signs of the LEs remain valid invariant indicators of chaos. Since we have extended the 
use of the LEs as a valid measure of chaos to include unbounded orbits, this conclusion is general: 
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it applies to both inertial and non-inertial reference frames and does not involve the identification 
of privileged observers. These results account for properties inherent to relativistic observers, 
such as event horizon and spatial unboundedness, significantly extending our understanding of the 
relativistic invariance of LEs and chaos. 
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